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rhé parallélisme autorisé par la mécanique quantique permet d'effectuer des calculs d'une manière 
radicalement nouvelle. Un ordinateur quantique fondé sur ces principes pourrait résoudre certains 
^^oblèmes exponentiellement plus vite qu'un ordinateur classique. Nous discutons comment un 
çMdinateur quantique réaliste peut simuler une dynamique complexe, en particulier les systèmes 
(réactiques où les fautes de l'ordinateur classique croissent exponentiellement vite. 



Prologue 
O 

"n se servit un verre de vin 
me, tira les rideaux et s'al- 
hgea pour réfléchir. Les équa- 
tions de la théorie du chaos ne 
jmsaient aucune référence au mi- 
^u physique dans lequel se déploy- 
&ent leurs manifestations; cette 
^Jquité leur permettait de trouver 
Etes applications en hydrodynamique 
l^pime en génétique des popula- 
^^^ns, en météorologie comme en 
S^iologie des groupes. Leur pou- 
voir de modélisation morphologique 
était bon, mais leurs capacités pré- 
dictives quasi nulles. A l'opposé, 
les équations de la mécanique quan- 
tique permettaient de prévoir le 
comportement des systèmes micro- 
physiques avec une précision ex- 
cellente, et même avec une préci- 
sion totales si l'on renonçait à tout 
espoir de retour vers une ontolo- 
gie matérielle. Il était au moins 
prématuré, et peut-être impossible, 
d'établir une jonction mathéma- 
tique entre ces deux théories." (M. 
Houellebecq, Les particules élémen- 
taires, Flammarion, 1998). 

INTRODUCTION 

Le problême soulevé ci-dessus 
devient encore plus frappant si on 



imagine que la dynamique chao- 
tique est simulée par un ordina- 
teur quantique. En effet, depuis une 
quinzaine d'années, la possibilité 
d'utiliser la mécanique quantique 
pour effectuer des calculs est de 
plus en plus discutée dans la com- 
munauté scientifique. Dans les an- 
nées 1980, Richard Feynman a pro- 
posé d'utiliser des éléments quan- 
tique pour simuler de manière effi- 
cace le comportement de systèmes 
quantiques. En effet, un système 
quantique à n corps évolue dans un 
espace de Hilbert dont la dimen- 
sion croît exponentiellement avec 
n. Par exemple, un système de n 
spins pouvant prendre deux valeurs 
évolue dans un espace de dimension 
2". Ceci oblige un ordinateur clas- 
sique à effectuer un nombre énorme 
d'opérations pour simuler le com- 
portement de ces systèmes même 
pour des n modérés. Par contre si 
les opérations sont physiquement 
effectuées par des éléments quan- 
tiques, il est possible de faire le 
même calcul avec des ressources 
bien plus petites. La question se po- 
sait donc de construire une théorie 
de l'information et une algorith- 
mique en tenant compte de ces 
effets quantiques. Ceci a été ré- 
solu ces dernières années, et on 



sait maintenant décrire un modèle 
théorique d'information quantique, 
fondée sur la notion de "qubit- 
s" (c'est-à-dire des systèmes quan- 
tiques à deux états |0) et |1)) sur 
lesquels on agit par des transforma- 
tions unitaires conservant la prob- 
abilité. Le qubit est l'équivalent 
quantique du bit, et décrit l'unité 
d'information quantique de la même 
manière que le bit en information 
classique. Néanmoins, il représente 
un objet essentiellement différent, 
pouvant varier sur des ensembles 
continus de valeurs, mais n'en four- 
nissant qu'une en cas de mesure. 
De plus, plusieurs qubits peuvent 
être intriqués, c'est à dire présen- 
ter des corrélations entre eux ab- 
sentes en information classique. Un 
exemple bien connu d'effets de l'in- 
trication apparaît dans le paradoxe 
Einstein-Podolsky-Rosen. En effet, 
pour deux qubits dans l'état in- 
triqué (|00) + |11))/V2, la mesure 
d'un qubit influera l'état de l'autre 
quelle que soit la distance entre eux. 

Un ordinateur quantique est 
généralement vu comme un ensem- 
ble de qubits sur lesquels on agit par 
des transformations unitaires. Ces 
transformations s'écrivent comme 
produits de matrices unitaires élé- 
mentaires appelées portes quan- 
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tiques et agissant sur un ou deux 
qubits (voir encadre 1). Le principe 
de superposition permet alors de 
créer un état quantique regroupant 
un nombre exponentiel d'états com- 
putationnels et d'agir sur eux tous 
en une seule opération. Par ex- 
emple, pour n qubits initialement 
dans l'état |00...0), l'application de 
n portes d'Hadamard (voir en- 
cadré 1) permet de créer l'état 
2-n/2 \x) qui regroupe tous 

les états computationnels de |00...0) 
à |ll...l). Des travaux ont permis 
d'établir que cette "parallélisation" 
massive due au nombre exponen- 
tiel d'états multi-qubits dans l'es- 
pace de Hilbert pouvait permet- 
tre d'accélérer de manière spectac- 
ulaire la vitesse de résolution de 
certains problèmes, d'abord dans la 
simulation de la mécanique quan- 
tique, mais aussi en d'autres cas 
comme la décomposition en facteurs 
premiers d'un grand nombre (ceci 
ayant des applications fondamen- 
tales en cryptographie). Une très 
grande activité se développe donc 
dans le monde en ce moment, pour 
essayer de construire un tel ordi- 
nateur, un problème technologique- 
ment très ardu où plusieurs voies 
sont exploitées pour la réalisation 
physique de qubits (voir l'article de 
M. Brune et J. M. Raimond dans ce 
numéro). À présent des algorithmes 
simples ont été implémentés dans 
des systèmes comportant jusqu'à 
sept qubits en utilisant les spins nu- 
cléaires de molécules par résonance 
magnétique nucléaire (RMN) ou des 
ions piégés. 

L'algorithme le plus célèbre 
est celui développé par Peter 
Shor qui permet de factoriser un 
grand nombre en facteurs pre- 
miers exponentiellement plus rapi- 
dement que n'importe quel al- 
gorithme classique connu à ce 
jour. En effet, le meilleur algo- 
rithme classique nécessite un nom- 
bre d'opérations qui croit comme 
0(exp[(logM)i/3(ioglogM)2/3]) 
où M est le nombre à factoriser, 
alors que l'algorithme de Shor re- 
quière O(n^) opérations quantiques 



élémentaires où n est le nombre 
de bits de M (n w logaAf). Cet 
algorithme a en fait une grande im- 
portance pratique, car bien que la 
factorisation soit difficile, l'opéra- 
tion inverse demande seulement 
O(n^) opérations pour reconstru- 
ire le nombre une fois les facteurs 
connus. Cette propriété d'asymétrie 
en fait la base du protocole cryp- 
tographique RSA utilisé partout 
dans le monde. En effet, l'opération 
facile permet de crypter aisément, 
et l'opération inverse difficile pro- 
tège le secret du message. Un ordi- 
nateur quantique permettrait donc 
de briser la plupart des codes de 
ce type. Néanmoins, les nombres 
utilisés par ces codes sont déjà très 
grands, et un ordinateur quantique 
d'au moins 1000 qubits serait néces- 
saire pour les factoriser suivant la 
méthode de Shor. Etant donné la 
taille des processeurs quantiques en 
construction, il est donc important 
de trouver des algorithmes quan- 
tiques dont l'utilité se manifesterait 
pour des nombres de qubits moins 
importants, par exemple quelques 
dizaines. 

Un exemple de ce type d'algo- 
rithme est représenté par la simu- 
lation de systèmes chaotiques, clas- 
siques ou quantiques. En effet, 
les systèmes classiques chaotiques 
présentent des propriétés d'instabil- 
ité exponentielle locale. Ceci im- 
plique que deux trajectoires ini- 
tialement voisines s'éloignent expo- 
nentiellement vite l'une de l'autre, 
ce qui rend les simulations de tels 
systèmes très difficiles. Pour cette 
raison, de tels systèmes sont de 
bons candidats de problèmes pour 
lesquels un ordinateur quantique 
pourrait être utile. Il est donc im- 
portant d'étudier la simulation de 
tels systèmes sur des ordinateurs 
quantiques. Un point crucial est 
également de comprendre comment 
les erreurs dues aux imprécisions 
et imperfections expérimentales af- 
fectent la précision du calcul quan- 
tique, dans le régime d'extrême in- 
stabilité représenté par le chaos 
classique (cf prologue) . 



ALGORITHMES QUANTIQUES 

Un algorithme classique est une 
suite d'instructions que l'on donne 
à l'ordinateur et qu'il doit ef- 
fectuer successivement en manipu- 
lant les bits de sa mémoire (qui va- 
lent ou 1). Un algorithme quan- 
tique agit de manière similaire sur 
des qubits, i.e. des systèmes quan- 
tiques à deux états |0) et |1). Si 
n tels qubits sont mis ensemble, le 
système total à tout instant sera 
décrit par une fonction d'onde, qui 
est un vecteur dans un espace de 
Hilbert de dimension 2". L'évolu- 
tion de ce système correspond à des 
transformations unitaires dans cet 
espace. Toute transformation uni- 
taire est licite, mais pour écrire 
un algorithme quantique, on re- 
quière d'écrire les transformations 
du système en fonction de transfor- 
mations élémentaires appelées aussi 
portes (voir encadré 1). Un algo- 
rithme quantique spécifie donc un 
état initial facile à préparer, puis 
une suite de transformations élé- 
mentaires (portes) à appliquer à cet 
état, et précise ensuite quelle infor- 
mation retirer de la mesure effec- 
tuée sur l'état final du système. 

L'algorithme de Shor dont il a 
été parlé précédemment suit cette 
structure générale. Il utilise le fait 
que le problème de factorisation 
peut être ramené au problème de 
trouver la période d'une certaine 
fonction /. Les qubits de l'ordina- 
teur quantique sont regroupés en 
registres, l'un d'entre eux contenant 
toutes les valeurs de a; de à A'' — 
1. Si ce registre contient n qubits, 
N — 2" et X peut ainsi varier 
de à 2" — 1. TV est choisi en- 
tre et 2M2, où m est le nom- 
bre à factoriser. Shor a alors mon- 
tré qu'il était possible de calculer 
en parallèle toutes les valeurs de 
f{x) de manière efficace sur un 
deuxième registre, permettant de 
passer de l'état N-^^^ J2x=o \^)\^) 
à iV-V2^^-^i|^)|/(^)). g. 

mesure les qubits du deuxième reg- 
istre, on fixe la valeur de f{x) = u, 
et la nouvelle fonction d'onde est 



une somme sur toutes les valeurs 
de X ayant pour image u par /. 
Mesurer le premier registre à ce 
stade ne donnerait pas d'informa- 
tion utile. On applique donc une 
transformée de Fourier quantique 
(QFT) (voir encadré 1) à cet état, 
obtenant une fonction d'onde finale 
concentrée autour des multiples de 
N/T, où T est la période. La déter- 
mination d'un de ces nombres par 
la mesure des qubits permet de cal- 
culer T. Un élément très important 
de cet algorithme est l'utilisation 



de la QFT, qui peut être effectuée 
en O(n^) opérations sur un vecteur 
de taille 2". La complexité algo- 
rithmique de l'algorithme de Shor, 
mesurée par le nombre de portes 
nécessaires pour factoriser un nom- 
bre M comportant n bits (n ~ 
log2 M), est de l'ordre de . Ce ré- 
sultat est obtenu grâce à l'utiHsa- 
tion de plusieurs propriétés spéci- 
fiquement quantiques : la possibil- 
ité d'agir sur tous les états compu- 
tationnels en même temps (principe 
de superpositions) puis (par l'action 



de la QFT) de les faire interférer de 
manière constructive (interférences 
quantiques). Ceci est impossible sur 
un ordinateur classique, et explique 
pourquoi l'algorithme de Shor est 
plus efficace que n'importe quel al- 
gorithme classique connu. Il est pos- 
sible d'implémenter l'algorithme de 
Shor sur un ordinateur classique, 
mais dans ce cas le nombre d'opéra- 
tions classiques est exponentielle- 
ment grand et l'efficacité quantique 
disparaît. 



Encadré 1 

PORTES QUANTIQUES 



Un ordinateur quantique est un ensemble de n qubits 
sur lesquels on agit par des transformations unitaires 
bien choisies. Toute transformation unitaire de 
l'espace de Hilbert de dimension 2" peut s'écrire 
comme combinaisons de transformations locales, 
agissant sur seulement quelques qubits, appelées portes 
quantiques. Un ordinateur quantique réel devra 
permettre l'implémentation répétée de quelques-unes de 
ces transformations élémentaires, et un algorithme 
quantique doit fournir la suite de pwrtes nécessaires 
pour pai'venir à l'état final. Quelques exemples usuels : 

-porte d'Hadamard s'appliquant à un qubit 
|0)-(|0) + |l))/x^;|l)-(|0)-|l))/x/2; 
-porte de phase s'appliquant à un qubit |0) — > |0) ; 

-controUed not ou CNOT s'appliquant à deux 

qubits : |00) ^ |00) ; |01) ^ |01) ;|10) ^ |11) ; 

|11) ^ |10) ; le deuxième qubit est inversé si le premier 

est dans l'état |1) ; 

-controUed controUed not ou porte de Toffoli 

s'appliquant à trois qubits : le troisième qubit est 
inversé si les deux premiers sont dans l'état |1). 

Certains ensembles de portes sont suffisant à eux seuls 
pour construire toutes les transformations unitaires : 
par exemple CNOT combiné aux transformations à un 
qubit. On parle alors de portes universelles. 

Une partie d'un algorithme quantique s'écrit donc 
comme une suite de portes, dont le nombre quantifie la 
complexité quantique du processus. Cette notion ne 



dépend pas du choix des portes élémentaires, 
puisqu'elles s'écrivent toutes les unes en fonction des 
autres. Un processus est polynômial si le nombre 
d'opérations nécessaires est une puissance du nombre 
de bits/qubits manipulés, exponentiel s'il en est une 
exponentielle. On peut vérifier que les opérations 
arithmétiques apprises à l'école primaire sont toutes 
des algorithmes polynômiaux. L'ordinateur quantique 
peut de plus additionner ou multiplier plusieurs 
nombres en parallèle. 

Une transformation unitaire générale est donnée par 
une matrice NxN et nécessite 0{N) portes 
élémentaires. Cependant, certaines transformations 
impwrtantes se décomposent en un nombre polynômial 
de portes. Un exemple de telle transformation, utilisée 
dans beaucoup d'algorithmes quantiques est la 
transformée de Fourier quantique (QFT). Elle utilise 
n qubits pour transformer un vecteur de taille 2" par : 

ELô^afcl^) — ' EL(î'^(ELô^exp(27riA://2")afe)|0. 
Elle peut s 'écrire au moyen des transformations 
élémentaires Hj (porte d'Hadamard appliquée au qubit 
j) et Bjk (porte à deux qubits appliquée aux qubits j et 
k et caractérisée par |00) |00) .- 
|01) ^ |01) ;|10) ^ |10) ; |11) ^ exp(ï7r/2'=-J)|ll);. 
On peut vérifier que la séquence H"^]^ [(n^^^^-^i3jfc)iJj] 
effectue bien la transformation de Fourier d'un vecteur 
de taille 2" en n{n + l)/2 opérations. 

En pratique, le choix des portes universelles dépend de 
l'implémentation expérimentale. 



D'autres algorithmes ont été 
aussi développés ; en particulier l'al- 
gorithme de Grover permet de faire 
une recherche dans une base de 
donnée non ordonnée de tafile N , 



avec un gain quadratique {0{\fN) 
opérations au lieu de 0{N) clas- 
siquement). Des algorithmes ont 
également été construits permet- 
tant de simuler efficacement cer- 



tains systèmes quantiques, en ac- 
cord avec la suggestion originale de 
Feynman. 
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Simulation quantique du chaos 
classique 

Il peut paraître naturel de pou- 
voir simuler la mécanique quan- 
tique efficacement sur un ordinateur 
quantique, bien qu'en fait une telle 
tâche ne soit pas évidente a con- 
crétiser sous forme d'un algorithme 
pratique. Il est plus surprenant 
d'imaginer qu'un ordinateur quan- 
tique puisse simuler efficacement la 
mécanique classique. Pourtant, les 
problèmes présentant du chaos clas- 
sique sont très difficiles à simuler 
sur un ordinateur classique. En ef- 
fet, l'instabilité exponentielle con- 
duit à une croissance exponentielle 
au cours du temps de la moindre 
imprécision sur la distribution clas- 
sique initiale. Un exemple d'appli- 
cation chaotique très connu est l'ap- 
plication du chat d'Arnold : y = y + 
x{modL)]x = X + ?;(modl), où y, x 
désignent les variables après une 
itération. Le mouvement a lieu dans 
l'espace de phase {x, y) sur un tore 
de taille L (entier) dans la direction 
y et 1 dans la direction x. Des ré- 
sultats mathématiques ont permis 
de prouver que la dynamique de ce 
système possède les caractéristiques 
du chaos très développé, avec par 
exemple un exposant de Lyapounov 
positif /i « 1. Ceci caractérise une 
divergence rapide d'orbites voisines 
dans tout l'espace de phase, la dis- 
tance entre elles augmentant avec 
le temps comme exp(/ii). En rai- 
son de cette instabilité, les erreurs 
d'arrondi dues à la précision finie 
d'un ordinateur classique vont se 
propager très vite dans le système. 
Par exemple, pour le Pentium III en 
double précision l'orbite sera com- 
plètement modifiée après un nom- 
bre d'itérations t = 40. Même si 
la dynamique exacte est réversible, 
l'existence de ces fautes couplées 
à l'instabilité du chaos détruit la 
réversibilité de l'évolution du sys- 
tème dans le temps. 

Grâce au parallélisme quan- 
tique, un ordinateur quantique peut 
simuler l'évolution d'un nombre ex- 
ponentiel de trajectoires classiques 
en un temps polynômial. Ceci per- 



met de calculer de manière fiable 
des quantités globales du système 
en dépit de l'instabilité exponen- 
tielle qui amplifie très vite toute im- 
précision. Trois registres sont néces- 
saires, deux spécifiant les valeurs de 
X et y et un servant d'espace mé- 
moire temporaire. En effet, on peut 
construire un état initial propor- 
tionnel à > IVj > |0 > 
où les aij valent ou 1, 1 < 
hj < 2". Il est possible d'effectuer 
les deux additions conduisant à l'é- 
tat Y,aij\2x.i + y.j > \xt + yj > 
|0 > de manière parallèle en seule- 
ment 0{n) opérations, en utilisant 
uniquement n — 1 qubits dans le 
troisième registre pour les retenues. 
De cette façon, 0(2^") trajectoires 
classiques sont simulées en même 
temps par un ordinateur quan- 
tique comportant 3n — 1 qubits au 
total. D'autres applications chao- 
tiques peuvent être simulées de 
manière similaire en un nombre 
polynômial de portes quantiques. 
Même une dynamique dissipative 
menant à un attracteur étrange 
peut être simulée de manière effi- 
cace. Après t itérations, l'état quan- 
tique du système contient les posi- 
tions des itérés des 0(2^") points 
initiaux. Une mesure de tous les 
qubits à ce stade donnerait seule- 
ment un point de cette distribution, 
et ferait perdre l'efficacité quan- 
tique. Pour éviter cela, il est pos- 
sible d'appliquer une QFT et de 
mesurer par ce biais des propriétés 
globales d'un nombre exponentiel 
d'orbites, obtenant ainsi une infor- 
mation nouvelle qui n'est pas acces- 
sible efficacement classiquement. 

Un ordinateur quantique réaliste 
comportera nécessairement des im- 
précisions dues au couplage avec le 
monde extérieur (voir l'article de 
M. Brune et J. M. Raimond dans 
ce numéro) qui feront que les portes 
idéales seront remplacées par des 
portes approximées qui vont intro- 
duire des fiuctuations d'amplitude 
e dans les transformations unitaires 
associées. Il est important de com- 
prendre si la dynamique chaotique 
va entraîner une augmentation ex- 



ponentielle de ces fautes quantiques 
comme c'était le cas pour les fautes 
d'arrondi dans l'ordinateur clas- 
sique. La Fig. 1 montre un exemple 
de la dynamique donnée par l'appli- 
cation du chat d'Arnold simulée par 
l'ordinateur classique et par l'ordi- 
nateur quantique. 




Figure 1 - Dynamique de l'application 
du chat d'Arnold simulée sur un 
ordinateur classique (à gauche) et 
quantique (à droite), sur un réseau de 
2'^x2^ points avec une seule cellule 
(L = 1 dans l'application du chat 
d'Arnold). Première ligne : distribution 
initiale ; seconde ligne : distributions 
après t = 10 itérations ; troisième ligne : 
distributions à t = 20 , avec renversement 
du temps effectué à t = 10 ; dernière 
ligne : distributions à t = 400, avec 
renversement du temps effectué àt = 200. 
A gauche : une erreur classique de la 
taille d'une cellule (e = l/128j est 
effectuée au moment du renversement 
du temps uniquement; à droite : toutes 
les portes quantiques opèrent avec des 
erreurs quantiques d'amplitude e = 0.01 ; 
la couleur représente la probabilité 
de bleu (zéro) à rouge (maximale) ; n = 7 
ce qui fait 20 qubits au total. 

Dans le cas de l'ordinateur clas- 
sique, 10 itérations sont suffisantes 
pour que la faute minimale sur le 
dernier bit détruise la réversibilité. 
En revanche, dans le cas de l'or- 
dinateur quantique, la réversibilité 
est préservée avec une bonne pré- 



cision en présence de fautes quan- 
tiques avec une amplitude compa- 
rable. La précision du calcul quan- 
tique reste raisonnable pendant un 
nombre d'itérations tj ex l/(ne^). 
La raison physique de ce résultat ré- 
side dans le fait que chaque porte 
imparfaite transfère une probabil- 
ité de l'état exact vers d'autres 
états. Ce résultat souligne la nature 
très différente du comportement des 
fautes dans le calcul classique et 
le calcul quantique. La théorie des 
perturbations quantiques explique 
la stabilité par rapport aux erreurs 
quantiques. Par contre, les fautes 
classiques ne rentrent pas dans ce 
cadre car même une très petite faute 
classique est très grande du point de 
vue quantique (renversement d'un 
qubit), et se propage exponentielle- 
ment rapidement dans le calcul clas- 
sique comme dans le calcul quan- 
tique. Chaque ordinateur a donc son 
type naturel de fautes qui se com- 
portent très différemment dans le 
cadre de la dynamique chaotique. 

Ceci permet de donner un nou- 
vel éclairage au vieux problème de 
l'irréversibilité en mécanique statis- 
tique, objet de nombreuses contro- 
verses depuis celle entre Loschmidt 
et Boltzmann au XIXe siècle. En 
effet, cette irréversibilité se mani- 
feste malgré le fait que les équations 
de Newton décrivant le mouvement 
microscopique sont réversibles. Une 
légende veut que Loschmidt ait 
souligné ce fait à Boltzmann en lui 
demandant ce qui arriverait à sa 
théorie statistique si on inversait 
les vitesses de toutes les particules, 
pour de cette manière retourner 
depuis un état d'équilibre à un état 
initial hors d'équilibre. La réponse 
de Boltzmann aurait été : "alors es- 
saye de le faire". Ceci serait possible 
sur un ordinateur quantique pour 
un nombre macroscopique de par- 
ticules. En effet, on peut par ex- 
emple utiHser l'algorithme ci-dessus 
pour L ^ 1. Dans ce cas, le sys- 
tème est étendu dans la direction 
y et il est connu que le chaos in- 
duit une diffusion dans cette di- 
rection. L'évolution statistique de 
la distribution initiale de partic- 



ules est alors donnée par l'équa- 
tion de Fokker-Planck. Un exem- 
ple de l'évolution du deuxième mo- 
ment < > de la distribution est 
montré dans la Fig. 2. Même si la 
dynamique exacte est réversible, en 
présence de fautes classiques, même 
très petites, la diffusion statistique 
reprend rapidement après le ren- 
versement du temps, comme af- 
firmé par Boltzmann. En revanche, 
sur l'ordinateur quantique même 
imparfait l'évolution statistique re- 
vient vers l'état initial hors d'équili- 
bre. Ceci permettrait donc à l'ordi- 
nateur quantique d'inverser la flèche 
du temps pour des systèmes macro- 
scopiques avec un nombre de par- 
ticules pouvant atteindre le nom- 
bre d'Avogadro (6.022x10^3) avec 
seulement 125 qubits (L = 8). 
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Figure 2 - Croissance diffusive du 
second moment < y'^ > pour l'application 
du chai d'Arnold avec L = 8, simulée 
sur des ordinateurs classique (Pentium 
III) et quantique ("Quantium I"). Le 
fonctionnement de Quantium I est bien 
sûr lui-même simulé au moyen d'un 
ordinateur classique. A t = 35 toutes les 
vitesses sont inversées. Pour le Pentium 
III l'inversion est faite avec une précision 
e = 10^^ (ligne poinlillée) cl e = 
10^^ (ligne en tirets); 10*' orbites sont 
.simulées. Pour Quantium I, la simulation 
est effectuée avec 26 qubits (n = 7) (ligne 
pleine épaisse); chaque porte quantique 
opère avec du bruit d'amplitude c = 
0.01. La ligne droite montre la diffusion 
macroscopique théorique. 

Simulation quantique du chaos 
quantique 

Les résultats exposés ci-dessus 
concernent la simulation du chaos 
classique. Il est aussi important d'é- 
tudier la simulation de dynamiques 
quantiques complexes sur des or- 
dinateurs quantiques. Une classe 
de systèmes particulièrement in- 
téressante correspond aux applica- 
tions quantiques. Leur dynamique 



s'écrit de manière générale ip = 
Uip on ^ est la fonction d'onde 
et 'ip sa nouvelle valeur après ac- 
tion de l'opérateur d'évolution U = 
^-ihP/2^-KV{e)/n La variable 9 
représente une phase et £ = —id/dO 
correspond au moment angulaire. 
Le paramètre h représente la con- 
stante de Planck sans dimension. 
Pour de nombreux choix de po- 
tentiel V, la dynamique du sys- 
tème classique correspondant de- 
vient chaotique quand K augmente. 
Le cas V{6) = coaO correspond à 
l'appHcation standard de Chirikov 
ï = 1 + KsinO; 6 = + ï, où les 
barres correspondent aux nouvelles 
valeurs de (/, 0) après une itération 
(l'action classique est I = M). Cette 
application apparaît dans la mod- 
élisation de plusieurs phénomènes 
physiques comme le confinement 
d'un plasma, les trajectoires des 
comètes ou le mouvement de par- 
ticules chargées dans des accéléra- 
teurs. Pour K > 0, ce système 
intégrable k K = montre une 
transition vers le chaos qui suit 
le théorème de Kolmogorov-Arnold- 
Moser. Pour K > 1, la dynamique 
classique devient chaotique avec in- 
stabilité exponentielle des orbites et 
apparition d'une diffusion dans la 
direction du moment £ (similaire à 
la diffusion montrée dans la Fig. 2). 
Le système quantique associé, ap- 
pelé rotateur puisé, représente un 
modèle très important dans l'é- 
tude du chaos quantique. Il mod- 
élise aussi la physique des atomes 
de Rydberg dans un champ micro- 
onde, permettant des comparaisons 
avec les résultats expérimentaux. 
De plus, il présente un phénomène 
de localisation dû aux interférences 
quantiques, par lequel un paquet 
d'onde quantique initialement local- 
isé dans l'espace des moments verra 
son étalement stoppé alors qu'une 
distribution classique de particules 
diffuserait dans tout le système. Ce 
phénomène est similaire à la locali- 
sation d'Anderson dans les solides, 
si bien que son étude permet de 
comprendre également la physique 
des électrons dans un système dé- 
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sordonné. Ce modèle a été réalisé 
expérimentalement par l'équipe de 
Mark Raizen à Austin (USA) avec 
des atomes froids. Le cas V{9) = 
{9"^ — a?Y correspond à un poten- 
tiel de double puits chaotique, et 
permet d'étudier le phénomène d'- 
effet tunnel quantique assisté par le 
chaos. En effet, la mécanique quan- 
tique permet à un paquet d'onde 
localisé initialement dans un puits 
(6 — —a ou — a) de passer d'un 
puits à l'autre périodiquement. Ce 
phénomène est similaire aux oscil- 
lations du chat de Schrôdinger. La 
présence de chaos quantique ren- 
force cet effet tunnel de manière 
importante, comme il a été vu ex- 
périmentalement avec des atomes 
froids. 




Figure 3 - Evolution du chat de 
Schrôdinger (paquet d'onde) animé par 
l'ordinateur quantique pour l'application 
quantique avec potentiel de double puits. 
La probabilité de présence en 9 (axe 
horizontal —tt < 6 < tt) est représentée 
en fonction du temps t (axe vertical), de 
t = (en haut) à t = 180 (en bas). Les 
paramètres sont ici K = 0.04, a = 1.6, 
n = iir/N avec N = 2("-i). Le calcul 
quantique est effectué avec n = 6 qubits, 
les portes sont imparfaites avec une 
amplitude de bruit e = 0.02. La couleur 
représente la densité de probabilité, de 
bleu (minimale) à rouge (maximale). 



Sur un ordinateur quantique, 
la dynamique de ces systèmes 
peut être simulée avec un nom- 
bre polynômial de portes pour un 
vecteur de dimension exponentielle- 
ment large TV = 2". Une itération 
nécessite 0{n^) opérations pour le 
rotateur puisé et 0{n'^) pour le 
potentiel de double puits, contre 
O(A^logA^) opérations sur un ordi- 
nateur classique. L'élément essentiel 
de l'algorithme quantique est l'util- 
isation de la QFT qui permet de 
passer efficacement de la représen- 
tation 9 à celle en £. 

La Fig. 3 montre un exem- 
ple d'une oscillation d'un "chat de 
Schrôdinger" simulée par l'ordina- 
teur quantique en présence de bruit 
dans les portes. Ce bruit induit une 
décroissance des oscillations suivant 
une loi en e~^* avec F ~ e^n^ où e 
est l'amplitude du bruit. La décrois- 
sance à taux F peut être considérée 
comme une manifestation de la dé- 
cohérence, c'est-à-dire une perte de 
la cohérence quantique due au cou- 
plage avec l'extérieur du système. 
La dépendance polynômiale de F 
dans le nombre de qubits montre 
qu'il est possible d'effectuer des sim- 
ulations de systèmes quantiques de 
très grande taille sur des ordina- 
teurs quantiques réalistes. 

IMPERFECTIONS DE 
L'ORDINATEUR ISOLÉ 

Le bruit dans les portes consid- 
éré ci-dessus correspond à une dé- 
cohérence due à un couplage avec le 
monde extérieur. Cependant, même 
en l'absence d'un tel couplage un 
ordinateur quantique isolé contient 
des imperfections statiques. En ef- 
fet, la distance en énergie entre les 
deux niveaux peut varier d'un qubit 
à l'autre, dans un intervalle 6. De 
plus, un couplage résiduel d'ampli- 
tude J entre les qubits est tou- 
jours présent. En effet, une inter- 
action entre qubits est nécessaire 
pour réaliser les portes quantiques 
(voir encadré 1). A première vue 
il peut sembler que J est toujours 
plus grand que la distance A„ entre 
les niveaux voisins de l'ordinateur 



quantique isolé. En effet, la densité 
de niveaux augmente exponentielle- 
ment avec n dans un tel système 
quantique à n corps, donnant la loi 
A„ ~ Sx2~" où n est le nombre de 
qubits dans l'ordinateur. On pour- 
rait penser que pour J > A„ les 
niveaux sont mélangés par l'inter- 
action entre les qubits et les états 
propres de l'ordinateur sont forte- 
ment modifiés par rapport à ceux de 
l'ordinateur parfait. Il est clair que 
limiter les interactions résiduelles 
à de telles valeurs serait impossi- 
ble même pour quelques dizaines de 
qubits. Heureusement, le mélange 
de niveaux n'apparaît en fait que 
pour des valeurs de J beaucoup plus 
grandes, qui dépassent le "seuil de 
chaos quantique" J > S/n. Cette 
valeur est beaucoup plus grande 
que A„ (exponentiellement) , et est 
liée au fait que l'interaction est 
à deux corps et couple un niveau 
à autres niveaux au plus. Au 
dessus de ce seuil, les caractéris- 
tiques du chaos quantique apparais- 
sent, un nombre exponentiel d'états 
sont mélangés, les fonctions propres 
de l'ordinateur isolé deviennent er- 
godiques, et les niveaux d'énergie 
obéissent à une statistique de ma- 
trices aléatoires. 




Figure 4 - Fusion de l'ordinateur 
quantique induite par le couplage entre les 
qubits. La couleur reflète l'entropie des 
états propres qui caractérise le mélange 
des états parfaits, du bleu (entropie nulle) 
au rouge (entropie maximale). Sur l'axe 
horizontal est représentée l'énergie des 
états propres de l'ordinateur comptée 
depuis l'état de plus basse énergie. L'axe 
vertical correspond de bas en haut à des 
valeurs croissantes du couplage résiduel 
J/5, de à 0.5. L'ordinateur quantique 
contient ici 12 qubits. Le seuil de chaos 
correspond à J/<5 ~ 0.3. 



Un exemple de "fusion" de l'or- 
dinateur induit par les couplages 
résiduels est montré sur la Fig. 4. Ce 
processus commence au centre de la 
bande d'énergie où la densité d'états 
est maximale et touche peu à peu 
tout l'ordinateur quand le couplage 
J augmente. Il est clair qu'il est 
préférable d'opérer l'ordinateur en- 
dessous du seuil de chaos. Au-dessus 
du seuil, le chaos se développe sur 
une échelle de temps r ~ nJ/J^. 
Au delà de r, l'usage de codes 
correcteurs d'erreur quantiques sera 
nécessaire pour continuer le calcul. 
De tels codes ont été développés, 
permettant d'effectuer des calculs 
pour des temps très longs en cor- 
rigeant au fur et à mesure les er- 
reurs, pourvu que bruit et impré- 
cisions quantiques restent suffisam- 
ment petits. 



CONCLUSION 

Les résultats présentés montrent 
que des systèmes aussi complexes 
que ceux présentant du chaos peu- 
vent être simulés avec une bonne 
précision sur des ordinateurs quan- 
tiques réalistes, même avec un nom- 
bre modéré de qubits, et de manière 
beaucoup plus efficace que sur un 
ordinateur classique. A présent, il 
est possible de réaliser des ordina- 
teurs quantiques de jusqu'à 7 qubits 
avec la technique de RMN. En 
principe ceci permettrait de simuler 
l'évolution présentée sur la Fig. 3, 
mais le temps de décohérence doit 
encore être amélioré pour permet- 
tre d'appliquer un plus grand nom- 
bre de portes élémentaires. 
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